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[Fe/H] relations for c-type RR Lyrae Variables based upon 
Fourier Coefficients 
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ABSTRACT 

[Fe/H] - - P relations are found for c-typ e RR Lyrae stars i n glob ular clusters. 
The relations are analogous to that found by lJurcsik fc Kovacsl (|1996[ ) for field ab- 
type RR Lyrae stars, where a longer period correlates with lower metallicity values 
for similar values of the Fourier coefficient cj>^\. The relations obtained here are used 
to determine the metallicity of field c-type RR Lyrae stars, those within to Cen, the 
LMC and toward the galactic bulge. The results are found to compare favorably to 
metallicity values obtained elsewhere. 

Key words: stars: abundances - stars: variables: other - globular clusters: general 



1 INTRODUCTION 

For nearly one hundred years, Fourier functions have been 
used to determine the pu lsation period, P, of regular vari- 
ables. [Snnon fc Led (|l98lh were able to show that other use- 
ful information apart from the period could be derived from 
Fourier functions. Their format of the relation is commonly 
used: 



V(t) 



Ao + ^ Ai cos(iut + i 



(1) 



where the terms Ai and (f>i are the Fourier coefficients of 
the fit of degree n, and u = 2n/P. Generally the values 
of Ai and (f>i were ignore d when this function was used to 
determine the value of P. Simon fc Led fll98lh showed that 
these terms when combined in the following manner - Rij = 
Ai/Aj and <j)ij — j't'i ~ i'Pii could provide more information 
about pulsating variables than at first thought. Currently, 
the coefficients derived from the light curves of Cepheids 
and RR Lyrae have been used to provide information about 
pulsation modes and resonance effects, as well as physical 
characteristics of the stars such as mass, luminosity, and 
metallicity. 

Jurcsik & Kovacs (f996, hereafter JK96) surveyed field 
ab-type RR Lyrae (RRab) and derived a relationship be- 
tween P, 03i and [Fe/H] for these variables which has been 
used in a variety of applications, in particular to determine 
the metallicity of stars found in large scale surveys (Mor- 
gan, Simet & Bargenquast f998, hereafter MSB98), in glob- 
ular clust ers (Cacciari, Corwin fc Carney 2005) an d in other 
galaxies (|Di Fabrizio et al.ll2005r i. ISandagd (2004) describes 
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this relationship as an aspect of the Oosterhoff-Arp-Preston 
period-metallicity effect, in which the apparent shift of the 
value of P or 03i is dependent upon the influence metallic- 
ity ha s on horizontal branch morphology. iKovacs fc Walker! 
(2001) expanded upon the work of JK96 by investigating 
possible relationships between the physical parameters of 
RRab stars and other Fourier coefficients. They found sev- 
eral trends in the value of photometric colours and magni- 
tudes with the coefficients. 

The utilization of the Fourier coefficients of the c-type 
RR Lyrae variables (RRc) has no t been as extensive as that 
of the RRab stars. ISimonl (|f 989l ) used hydrodynamic mod- 
els of RRc stars to show relationships exist between the 
values of helium abundance (Y), luminosity, mass, P and 
031. In particular, he found that 03i is directly rel ated to 
the lu minosity-mass ratio, L/M 1 ' 81 , for these stars. [Simonl 
(1990) showed that the distribution of 03i with period would 
vary according to Oosterhoff group, with Oosterhoff I clus- 
ter stars having larger values of 03i for a given value of P 
than stars from Oosterhoff II clusters. Clement, Jankulak & 
Simon (1992) found strong evidence for the trend of 03i in- 
creasing with period in their observations of several globular 
clusters, and also noted the trend for lower metallicity clus- 
ters to have their 03i values shifted to longer periods than 
the higher metallicity clusters on a 03i-P diagram. This re- 
sult was also observed in the RRc variables in NGC 4590 
(M68) by Clement, Ferance & Simon (1993). The general 
nature of this aspect of the variation of 03i with metallicity 
has previously been used to estimate the metallicity of RRc 
stars in the OGLE s urveys of the galactic bulge (MSB98) 
and toward 47 Tuc ((Morgan fc DickersorJ 12000). however 
these results were at best only general approximations of 
[Fe/H]. 
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The lack of a metallicity relationship for the RRc stars 
analogous to that derived by JK96 is likely due to the 
scarcity of accurate metallicity and light curve data for field 
RRc variables. Generally metallici ties for field RR c stars are 
estimated using the AS method of |Prestonll|l959ll . with only 
a few metallicities derived from high resolution spectroscopy 



Table 1. Globular Clusters containing RRc Stars with reliable 
Fourier Coefficients. 



llButler et all 1 1982 



1 19971 ; ISolano et al 



.Lambert et al.l 1 19961 : iFernlev fc Barnea 
1 19971 ). In several cases, AS values ob- 



tained for a single RRc star can vary significantly, making 
the accuracy of the metallicity for such stars suspect. 

Fortunately, there is an abundance of data for globular 
cluster RRc variables, including values for cluster metallic- 
ity and a significant number of Fourier coefficients. We will 
use these data to derive relationships for RRc stars relating 
the Fourier coefficient <j)3i, pulsation period and [Fe/H]. The 
resulting relations will also be applied to several test cases, 
including field RRc stars, as well as variables in ui Cen, the 
Large Magellanic Cloud and toward the galactic bulge. 



2 METHOD 

The Fourier coefficients for globular cluster RRc variables 
were obtained from the Fourier coefficient website l|Morganl 
l2003h as well as several recent publications. The clusters 
used here are listed in Table [1] In order to avoid possible 
errors with conversion from one photometric system to an- 
other, only Fourier coefficients derived from V magnitudes 
were considered. Where ever possible the light curve data for 
the individual stars were examined and Fourier fits that were 
based upon sparsely sampled light curves or with poorly 
defined maxima/minima were excluded. Fourier coefficients 
with large uncertainties were also excluded. The sources for 
the Fourier coefficients, and the number of stars from each 
cluster ultimately used in this study are included in Ta- 
ble Q] as well. There are several sources for metallicity that 
could be used for globular clusters. Two metallicity scales 
that are frequently cited in the literature are those of Zinn & 
West (1984, hereafter ZW84) and Carretta & Gratton (1997, 
hereafter CG97). Both of these metallicity scales will be used 
here. For some clusters, the updated values from lZinnl (|l985l ) 
are used in place of the ZW84 values where available. When 
values were not available for clusters based on the system 
of CG97, the values found by Z W84 were converted to the 
CG97 scale using the relation of ICarretta et al.1 i|200ll ). The 
metallicity values and their corresponding uncertainties are 
also given in Table [T] 

Overtone pulsators such as RRc stars generally have a 
lower order fit (small value of n in equation [1} used for the 
light curve, which limits the number of Fourier coefficients 
available. As was out lined in the introduction and shown in 
IClement et~aH i| 19921 ) and MSB98, the coefficient that ap- 
pears to depend strongest upon [Fe/H] is the cj>zi term. This 
tren d is also observed in RRab stars as was noted by JK96 
and lClement fc Sheltonlll 19991) . Our study concentrated ex- 
clusively on this coefficient. The (j>z\ values for the 106 stars 
in our sample are plotted relative to P in Figure [T] with the 
coding in the diagram based upon the [Fe/H] values from 
Tabled] 

The relation between [Fe/H] - P - (f>3i found by JK96 for 
RRab stars is linear as is a simila r relation found bv lSandagei 
l|2004l ). It should be noted that ISandaed (|2004l ) used logP 



Cluster 


Symbol 


Source 


N 


[Fe/H] zw 


[Fe/H] CG 


NGC 6171 


X 


1 


7 


-0.99 ±0.06 


-0.97 ± 0.04* 


NGC 6362 


X 


2 


13 


-1.08 ±0.09 


-0.96 ± 0.01 


NGC 1851 


o 


3 


4 


-1.33 ±0.09 


-1.18 ± 0.05* 


NGC 5904 


o 


4 


14 


-1.40 ±0.06 


-1.11 ± 0.11 


NGC 6934 


A 


5 


6 


-1.54 ±0.09 


-1.32 ± 0.07* 


NGC 7089 


A 


6 


3 


-1.62 ±0.07 


-1.38 ± 0.06* 


NGC 5272 


A 


7 


17 


-1.66 ±0.06 


-1.34 ±0.06 


NGC 6333 


+ 


8 


5 


-1.78 ±0.15 


-1.52 ± 0.16* 


NGC 4147 


+ 


9 


8 


-1.80 ±0.26 


-1.55 ± 0.28* 


NGC 6809 


+ 


10 


5 


-1.82 ±0.15 


-1.57 ±0.17* 


NGC 4590 


V 


11 


14 


-2.09 ±0.11 


-1.99 ± 0.10 


NGC 7078 


V 


12,13 


10 


-2.15 ±0.08 


-2.12 ± 0.01 



Sourc es: 1 
ll200ll) . 3 



- IClement fc Sheltonl JlQQTt). 2 - lOIech et al.l 
IWalkeJ dl998h, 4 - IKaluznv et al] d2000l)~ 5 



iKaluznv. 01ech~fc**StanelJ d200ll), 6 - lAlcock et al l d2004l) , 
7 - ICacciari et al.l d2005f>7 8 - IClement fc Sheltonl dl999l). 9 
- IStetson, Catelan fc Smithl (|2005|), 10 - lOIech et alj lll999l) , 
11 - IWalkeJ dl994h. 12 - IClement fc Sheltonl jl996h . 13 - 
lArellano*l^rro. Garcia Lugo, fcl^ose nzweigf (2006) 
* - calculated using ICarretfa^t**alT^|2*0*0*lf i 
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Figure 1. Fourier coefficient 1/131 plotted versus period for the 
RRc stars in the globular clusters listed in Table [TJ The various 
symbols are based upon the cluster's value for [Fe/H]. 



rather than P with virtually the same degree of accuracy 
as that obtained with the JK96 relation. We found the use 
of a P term to be slightly better than the use of a logP 
term when comparing the residuals of the derived relations 
to the observed values. The form of the relationship between 
[Fe/H], P and (j>3i appears to depend upon the metallicity 
scale that is used, ZW84 or CG97. This is seen when the 
P values for the RRc stars in each cluster are increased by 
the same amount to place the </>3i values along a common 
distribution, as is illustrated in Figure [2] The dependence of 
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Figure 2. Fourier coefficient 03i plotted versus the shifted value 
of P. All variables within a given globular cluster are shifted by 
the same amount in P. 

the 03i values on P is apparent, while the amount that was 
added to the period of each cluster's RRc stars indicates the 
dependence of P on the cluster's metallicity. The amounts 
added to the periods of the RRc in each cluster in order to 
produce the distribution in Figure [2] are shown in Figure [3] 
For the CG97 scale, a linear relation appears to exist be- 
tween the shift in P and [Fe/H], while for the ZW84 scale a 
non-linear relation appears to be appropriate. In the case of 
the ZW84 data, the standard deviation of the points from a 
linear fit is cr = 0.127, while a quadratic fit has a standard 
deviation of a — 0.092. Applying a quadratic fit to the CG97 
data results in a function that is nearly identical to the lin- 
ear one, indicating that no improvement to the fit is made 
with a higher order function. Even though the CG97 scale 
does not appear to require a high-order non-linear relation- 
ship between [Fe/H] and P, those relations were nonetheless 
examined. 

In order to determine the best relationship between the 
three variables, all possible permutations were made be- 
tween them with functions of the form 

[Fe/H] = a P 2 + b P + c 0^ + d 3 i + e P0 3 i + /, (2) 

tested for quality of fit to the data. All possible combina- 
tions of zero and non-zero values for the coefficients a - f 
were examined using both the ZW84 and CG97 scales. Only 
relations which included at least one P and 03i term were 
considered. A summary of the resulting fits to the various 
formulae are shown in Table [2] This table shows the range of 
values for the sample standard deviation of the best fitting 
formulae to the 106 data points. The best solution for the 
ZW84 scale is one with all six terms used in equation(2] while 
the quality of the best solutions for the CG97 scale does not 
vary significantly when fewer terms are used in equation [5] 
This result is expected given the linear relation of the [Fe/H] 
values on P as is shown in Figure [3] for the CG97 system. 



Figure 3. The amount of shift in period used to produce the 
distribution shown in Figure [2] for each cluster plotted versus the 
cluster's value of [Fe/H] for the ZW84 (top) and CG97 (bottom) 
systems. The best fit quadratic and linear relations are drawn 
through the ZW84 data, while only a linear relation is shown for 
the CG87 system. 



Table 2. Sample Standard Deviations (it) for solutions to equa- 
tion [2] 



Metallicity 


Terms 


Range of a 


ZW84 


6 


0.145 




5 


0.152 - 0.161 




1 


0.160 - 0.170 




3 


0.162 - 0.174 


CG97 


6 


0.142 




5 


0.142 - 0.144 




4 


0.142 - 0.154 




3 


0.143 - 0.169 



Two criteria were used for determining which formulae 
would be the most useful for calculating RRc metallicities. 
First was the quality of the fit to the data, which is summa- 
rized in Table [2] The second criterion was the simplicity of 
the formula. This criterion was only relevant for the CG97 
solutions, where the quality of the formula varied insignifi- 
cantly when the number of terms used in the solution were 
changed. The best fit formula for equation [2] to the clusters 
in Table Q] using the Z W84 scale is 

[Fe/H]zw = 52.466P 2 - 30.075P + 0.131031 

+0.982031 - 4.1980 3 iP + 2.424 (3) 

which has a sample standard deviation of 0.145 dex. The 
best formula based upon the CG97 scale is 

[Fe/H]cc = 0.0348031 + O.1960 3 i - 8.507P + 0.367 (4) 

which has a sample standard deviation of 0.142 dex. 84% 
and 87% of the [Fe/H] values based upon the above for- 



4 S. M. Morgan, J. N. Wahl and R. M. Wieckhorst 



2 



_ Equation 3 




ii 


- — : 










- /^$k 

^- — — 1 — H 


<A V / 

-2' 




: 


_ Equation 4 


1 1 


— i — i — i — — 


— 1 — ^ — i^- — — ^ 


: / /£■ 

- -1/ _ j 2s x - 


/v 

1 75/ 

' 1 / 







0.2 0.3 0.4 0.5 



Period 

Figure 4. Lines of constant [Fe/H], based upon the fits to the 
ZW84 (equation [3} and the CG97 (equation |4t metallicity scales 
are shown, along with the RRc data from Figurefl] [Fe/H] values 
for each line are given in the lower left corner of each graph. 

mulae are within 0.2 dex of the ZW84 and CG97 cluster 
metallicity values respectively, while 58% and 62% of the 
[Fe/H] values are within 0.1 dex. The average difference be- 
tween the [Fe/H] values from ZW84 and CG97 and those 
based upon the above relations is approximately 0.14 dex, 
indicating that these relations provide values for the metal- 
licity that are within the range of uncertainty found using a 
variety of other methods. 

Lines of constant value of [Fe/H] based upon equations[3] 
and [4] are shown in Figure |4j along with the original data 
from Figure [T] The non-linear nature of equation [3] limits 
the range of values where the function is defined in Fig- 
ure [4] The lines follow the general metallicity trends seen 
in the two metallicity systems. Average metallicities for all 
of the variables in each cluster listed in Table [1] based upon 
equations [3] and [4] are shown in Figure [5] Error bars in Fig- 
ure[5]for the average calculated metallicities are taken as the 
standard deviation of the [Fe/H] values for the stars in each 
cluster. The average cluster metallicity based upon equa- 
tion [3] for the 12 clusters shown in Figure [5] varies from the 
ZW84 values with a mean of 0.08 dex, while the values based 
on equation [4] vary from the CG97 values with of mean of 
0.07 dex. NGC 7089 (M2) is the most displaced cluster in 
Figure [5] with two of its three stars having relatively high 
calculated metallicities compared to the values of ZW84 and 
CG97. This results in an average metallicity for NGC 7089 
that is 0.22 and 0.20 dex greater than the ZW84 and CG97 
values respectively. 

3 TEST OF THE [FE/H] RELATIONS 

Equations [3] and [4] were applied to several test cases to ex- 
amine the quality of the relationships in determining [Fe/H] 
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Figure 5. Average cluster metallicities based upon ZW84 and 
equation [3] (top), and CG97 and equation [4] (bottom) are shown, 
along with a line of value unity. 

in other environments. Field RRc stars were examined first. 
The number of field RRc stars with [Fe/H] values and good 
quality V magnitude light curves that could be used is very 
small, comprised of on ly 15 stars. [Fe/H] values for the stars 
were taken either from lFernlev fc Barnes! l|l997h or were cal- 
culated using published AS valu es and the A5* - [Fe/H] re- 
lation of lFernlev fc Barnesl l| 19971 ). Some stars have only one 
measured value of AS, while others have widely divergent 
values. The comparison of the [Fe/H] values from the liter- 
ature and equations [3] and [4] is shown in Figure [6] The esti- 
mated uncertainty of the individual values of [Fe/H] from the 
literature varies with each source, but typical uncertainties 
are 0.2 dex or less. This uncertainty is similar to that found 
in the derivation of equations [3] and [4] (approximately 0.15 
and 0.14 dex respectively). These uncertainties are displayed 
in Figure [6] The sample standard deviation of the average 
metallicity values from the literature and those derived us- 
ing our formulae is 0.41 dex for equation [3] and 0.42 dex 
for equation [4] There are several stars that are well removed 
from the unity relation in Figure [6] These are TV Boo, ST 
CVn and V487 Sco. TV Boo ([Fe/H] = -2.44) has a metal- 
licity value from the literature outside of the range used to 
derive equations and 2] which may explain its divergent 
value. ST CVn (P = 0.329025 days) has a relatively large 
values of 03i (5.11) for its value of P, which accounts for 
its abnormally high calculated values of [Fe/H]. V487 Sco's 
metallicity ([Fe/H] = —1.89) is based upon a single value 
for A5. When these stars are excluded, the sample stan- 
dard deviations reduce to 0.26 dex (equation [3J) and 0.15 
dex (equation [3}. 

The RRc stars in uj Cen were examined next. Data 
from the variables in this cluster were collected from the 
Fourier coefficient website ( Morgan 2003) and an initial data 
set was selected based upon the method of observation, 
where recent, CCD based coefficients were favored over pho- 
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Figure 6. Values of [Fe/H] from the literature for field RRc stars 
compared to values based upon equations \3\ (top) and[4](bottom). 
TV Boo, V487 Sco, and ST CVn are indicated. The line indicates 
a ratio of unity. 



Figure 7. The range of metallicity values for RRc stars in ui Cen, 
based upon equations [3] (top) and [4] (bottom). 



togra phic observations. These in cluded Fourier coe f ficient s 
from (Morgan fc Dickersonl (|2000T ). and lOlech et all (|2003h . 
As with the globular cluster variables used to derive equa- 
tions [3] and 31 the RRc stars in ui Cen were examined and 
poor quality fits or sparse light curves were excluded. This 
reduced the number of variables to 67, and equations [3] and 
|3] were applied to these stars. The resulting distribution of 
metallicity is shown in Figure [7] The relatively wide spread 
in metallicity is not surprisin g given the broad range of 
metallicity observed in cu Cen ((Sollima et all 120051 ) . Values 
that we calculated for [Fe/H] range from —2.04 to —0.97, 
with a mean of —1.65 ± 0.24 using equation [3] and a range 
of -2.22 to -0.83, with a mean of -1.44 ± 0.25 for equa- 
tion!?] The largest concentration of values is near —1.77 and 
— 1.50 for equ ations [3] and [4] respectively. These peaks are 
near the value [Sollima et al. (2005) found for the dominant 
metallicity population of uo Cen, [Fe/H] ~ — 1.7 dex. 

There have been several studies that have found the 
metallicities of individual stars within ui Cen using a var i- 
ety of methods, such as Caby photometry (iRev et al.ll2000h . 
high-res olution spectroscopy l|Sollima et aL " r2006l ) and line 
indices (|Gratton et al.ll2004h . These methods cover a total 
of 54 stars that also have good quality light curves and de- 
rived Fourier coefficients. The individual stars are compared 
to the metallicities based on equations [3] and U in Figure [5] 
The values based upon the ZW84 relation (equation [3]) have 
a sample standard deviation from the published metallicity 
values of 0.25 dex, while the CG97 values (equation Q have 
a deviation of 0.23 dex. Errors for individual stars are not 
shown in Figure [8] due to the crowded graph, however most 
values for errors from the literature are approximated at 0.2 
dex. Errors from equations [3] and [4] are again taken to be 
0.15 and 0.14 dex respec tively. 

lAlcock et al.l (|2004T ) provided l/-band Fourier coeffi- 





Figure 8. Individual RRc stars in u> Cen with metallicities from 
the literature, compared to those based upon equations \3\ (top) 
and [4] (bottom). The line indicates a ratio of unity. The typical 
errors associated with the data points arc shown. 



cients for 682 RRc in the Large Magellanic Cloud. Equa- 
tions [3] and [4] were applied to all of these stars and re- 
spective metallicity averages of [Fe/H] = —1.61 ± 0.40 and 
[Fe/H] = —1.42 ± 0.37 were found. There were significant 
deviations from these averages, with some positive metallic- 
ity values calculated. When stars with large errors (> 0.5) 
in their values of 03i were removed from the sample, the re- 
sulting average metallicity changes slightly, to —1.66 ± 0.29 
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4 CONCLUSIONS 

Two formulae were derived that show the relationship be- 
tween the Fourier coefficient <^3i, pulsation period and 
[Fe/H] for RRc stars in 12 globular clusters. The for- 
mulae are based o n the widely used metall i city s cales of 
IZinn fe West] i| 19841 ) and lCarretta fe Grattonl (| 19971 ). These 
relations (equations [3] and were able to provide reliable 
estimates for the value of [Fe/H] for RRc stars in other envi- 
ronments, including uj Cen, the LMC, the galactic bulge and 
field RRc stars. Even though these results are encouraging, 
we do realize that the metallicity relations found here will no 
doubt evolve as more high quality light curves for globular 
cluster RRc stars are made available. In particular, values 
of (j>3i for RRc stars in clusters with very high or very low 
metallicities would be useful to refine and improve the for- 
mulae derived here, which were based upon cluster [Fe/H] 
values between approximately —1 and —2. At the present 
time it is hoped that the formulae presented here will be of 
use to others investigating the characteristics of RR Lyrae 
variables. 



Figure 9. Individual RRc star s in the LMC with metallicities 
taken from [Gratton et al. I ll2004l) , compared to those based upon 
equations [3] (top) and l4l (bottom). The line indicates a ratio of 
unity. 



(equation [3j and —1.45 ± 0.31 (equation |H) . The average 
[Fe/H ] values are very similar to those of iGratton et al.l 
( 2004), who found a value of -1.48 dex for 101 RR Lyrae, 
and also very similar to the met allicity of 23 LMC RR Lyrae 
found bv lBorissova et al] (120041 ) . [Fe/H] = -1.46. It is pos- 
sible to do a s t ar-by -star comparison to 14 LMC stars from 
IGratton et al] l|2004l ) , who calculated metallicities using line 
indices. The statistical comparison to the 14 stars' metallic- 
ities results in sample standard deviations of approximately 
0.28 for both equation [3] and [4j with the star-by-star com- 
parison illustrated in Figure [5] It should be noted that the 
met allicities of ind i vidual stars in the LMC and ui Cen from 
the IGratton et al I (120041) study were calibrated using the 
[Fe/H] scale of iHarris! l| 19961) , which is typically closer in 
value to the ZW84 system. 

The Fourier coefficients for 60 RRc stars observed in 
the direction of the galactic bulge were derived by MSB98. 
These values were obtained from J-band photometry, and 
were transformed to the I^-band using the relations from 
MSB98. The average metallicities obtained using equations[3] 
andg]were -1.03±0.50 and -0.98±0.37 respectively. There 
are several notable outliers in the sample though, some of 
which have been previous noted by MSB98 as likely having 
unusual metallicity values. These include BW9 V38, and 
BW11 V55, both of which have abnormally high values of 
[Fe/H], and BW1 Vll, BW2 V8, BW2 V10, BW4 V46, BW7 
V30, BW10 V45 and BW11 V34, all with at least one value 
of [Fe/H] < —1.75. When these 9 stars are excluded, the 
metallicity averages become —0.97 ± 0.36 (equation [3]) and 
—0.91 ± 0.22 (equation Q. These metallic i ties c ompare fa- 
vorably to the value obtained by ISmoled ((2005) for RRab 
stars in the galactic bulge ([Fe/H] = -1.04 ± 0.03). 
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